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ABSTRACT 


Finite  state  space  semi-Markov  processes  find  application  in  many  areas.  Often  in¬ 
terest  centers  on  whether  or  not  the  process  has  hit  a  particular  state  before  a  time  t. 
This  thesis  reports  results  of  a  simulation  study  of  the  small  sample  behavior  for  three 
estimators  of  the  survival  probability  of  a  first  passage  time  for  a  semi-Markov  process 
using  censored  data. 
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I.  INTRODUCTION 


Finite  state  space  semi-Markov  models  find  applications  in  a  variety  of  fields  such 
as  queueing  theory,  reliability,  and  clinical  trials.  Often  interest  in  the  application  of 
these  models  centers  on  the  distribution  of  the  first  passage  time  to  a  state  or  a  set  of 
states  representing,  for  example,  the  lifetime  of  a  system  or  the  end  of  a  busy  period  of 
a  server.  Suppose  that  the  observations  of  the  path  of  the  semi-Markov  process  are  all 
that  is  known  about  the  process.  The  problem  is  to  estimate  the  probability  that  the  first 
passage  time  has  not  occurred  before  time  t. 

Censored  data  problems  arise  frequently  in  medical,  and  also  in  engineering  system 
reliability  applications.  For  example,  in  medical  survivorship  studies  some  subjects  may¬ 
be  lost  to  follow-up,  or  available  data  may  be  analyzed  before  all  subjects  have  expired. 
In  the  equipment  reliability  context,  observed  units  may  still  be  in  operation,  perhaps 
after  several  previous  failures,  at  the  time  of  the  analysis. 

Three  possible  estimators  will  be  considered.  The  three  estimators  use  different 
amounts  of  information  concerning  the  process.  One  estimator  uses  only  the  observed 
first  passage  times.  Another  estimator  makes  parametric  assumptions  concerning  the 
sojourn  time  distribution  and  uses  maximum  likelihood.  A  third  approach  uses  an  ex¬ 
ponential  approximation  to  the  probability  and  empirical  distributions  to  estimate  the 
sojourn  time  distributions. 

The  three  estimators  were  investigated  for  a  specific  semi- .Markov  process  with 
uncensored  data  in  Kim[Ref  1]  and  Jacobs[Ref  2].  Results  of  a  simulation  study  of  the 
three  estimators  using  censored  data  are  reported  in  Gallagher  [Ref  3].  The  emphasis  in 
this  latter  study  is  on  the  behavior  of  the  point  estimates  with  mean  biases  and  standard 
errors  being  given. 

In  this  thesis  the  in  vestigation  of  the  behavio’-  of  the  three  estimators  with  censored 
data  is  continued.  The  emphasis  here  is  on  confidence  intervals  for  the  three  estimators. 
In  Chapter  2,  the  three  estimators  are  described  and  the  respective  confidence  interval 
procedures  considered  introduced.  Chapter  3  contains  the  details  of  the  simulation  ex¬ 
periment  and  its  results.  Finally,  conclusions  from  the  study  are  given  in  Chapter  4. 
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11.  NATURE  OF  PROBLEM 


A.  PROBLEM 

The  semi-Markov  process  model  considered  is  as  follows.  Suppose  we  observe  N 
individuals.  Let  A',(/)  be  the  state  of  the  /'*  individual  at  time  t.  We  will  assume 
{T,(/);r>:0}  /  =  1,2, 3,...., A.  are  independent  identicaUy  distributed  semi-Markov  proc¬ 
esses  having  the  same  probability  law  as  {X,\  r  ^  0}.  The  process  {X;,  r  ^  0}  is  a  semi- 
Markov  process  with  three  states  {0, 1,2}.  The  individuals  start  at  t  =  0  in  state  1. 
Upon  leaving  state  1,  the  process  transitions  to  state  0  with  probability  $  and  to  state  2 
with  probability  1  —  6.  From  state  2  the  process  transitions  to  state  I  with  probability 
I.  State  0  is  an  absorbing  state.  The  first  passage  time  to  state  0  will  be  referred  to  as 
the  time  of  death.  The  N  individuals  are  censored  independently.  The  censoring  times 
are  exponentially  distributed  with  a  mean  of  -p.  The  entire  path  of  transitions  and 
sojourn  times  are  observed  until  the  ume  of  censoring  or  death. 

Let 


and 


SW  =  P{D>i} 

where  D  is  the  time  of  death  (or  entrance  to  state  0).  The  problem  is  to  estimate  the 
survival  probability  P{D  >  /}  with  the  censored  data  of  N  individuals. 

B.  ESTIMATORS 

1.  Kaplan-Meier  Estimator 

A  non-parametric  estimator  of  the  distribution  function  for  censored  data  is  the 
Kaplan-Meier  estimator(K.M.E.)  which  is  often  called  the  product  limit  estimator  [Ref. 
4).  Let  (/,,  Uj,  ...,(/„  be  independent  identically  distributed  random  variables  with  distrib¬ 
ution  G  having  a  density  function.  Let  F,,  Fj, ...,  F„  be  independent  identically  distrib¬ 
uted  times  to  censoring  with  a  continuous  distribution  function.  Let 

Zi  =  min(f;,-,  Vj), 


and 


2 


if  U^<Vi 
Otherwise. 


The  Z,  are  the  observed  times  and  S,  is  an  indicator  of  whether  or  not  the  observation 
is  censored.  Let  Z(i)  ^  Z(2,  <  ...  <  Z|j„  be  the  order  statistics  of  {Z,}  and  be  the  corre¬ 
sponding  values  of  {<),}.  It  is  assumed  that  there  will  be  no  ties  since  the  underlying 
distribution  functions  are  continuous.  The  Kaplan-Meier  estimate  of  the  survival  func¬ 
tion  S{t)  =  {1  -  G(r)}  is 


S{i)  = 


flfrai 


<r| 

0 

Undefined 


if 

if  ^  ^  ^{n)  ^  ^(n)  ~  ^ 
if  ^  >  ^(n)  ^  ^{n)  —  1  • 


The  variance  of  S{i)  is  given  approximately  by 


VarlS{t)’\  ^  [5(/)]' 


y _ ^ 

{n~i){n-i- 


(':Z,o  <r] 


[Ref  4:  p.  464).  The  Kaplan-Meier  estimator  using  the  death  times  of  the  N  individuals 
will  be  denoted  by  Pk{D  >  /}. 

If  the  Kaplan-Meier  estimate  is  undefined,  we  investigate  the  effect  of  two 
methods,  defined  as  MOD  1  and  MOD  2,  to  make  the  Kaplan-Meier  estimate  honest. 
MOD  1  and  MOD  2  are  defined  as  follow.  In  MOD  1,  the  remaining  mass  of  the  esti¬ 
mated  survival  function  is  assigned  to  the  last  datum  Z^jfwhich  is  censored).  In  MOD 
2,  if  the  last  k  data  points  are  censored  the  remaining  mass  of  the  survival  function  is 
distributed  equally  among  the  k  data  points.  For  example,  if  the  estimated  survival 
probability  at  the  last  uncensored  point  is  0.2  and  there  are  additional  two  data  points 
w'hich  are  censored,  then  each  of  these  additional  points  is  assigned  a  mass  of  — r-  . 


2.  Maximum  Likelihood  Estimator 


In  this  subsection,  the  maximum  likelihood  estimator(M.L.E.)  will  be  given  for 
the  special  case  in  which  the  sojourn  time  in  state  /  is  exponentially  distributed  with 
mean -^(/=  1,2). 

Let  R,j  be  the  number  of  transitions  from  state  /  to  state  j  for  one  individual. 

The  log  likelihood  function  for  an  individual  is  ‘ 

/  =  ln(l  —  0)  +  /?io  In  d  +  In p2  "b  (^lo  ^12)  Pi  ~  Pi ^1  ~  Pi^i  (2-3) 


where  (/  =  1,2)  is  the  total  time  spent  in  state  /  before  death  or  censoring. 

The  maximum  likelihood  estimators  using  the  data  from  all  N  individuals  are 


MO 


^10  +  ^12 


Pi  = 


^10  "b  ^12 


Tx 


A  _  ^21 

7; 


A  U  death  times  that  are  censored 
c  = - r - r - 

r, +  r, 


(2.4) 


(2.5) 


(2.6) 

(2.7) 


where 


=  (2.8) 

and 

(2.9) 

with  /?,/«)  being  the  number  of  transitions  from  i  to  j  for  the  individual  and  Tin)  being 
the  total  time  spent  in  state  i  before  death  or  censoring  for  the  /'*  individual. 

To  obtain  asymptotic  variances  for  these  estimators,  note  that 
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(2.10) 


Further 


-dh  _  ^10  ,  ^12 

ee^  [i-ef 

—8^1  _  ^10  ^12  , 

2  “  2  ’ 

0P\  Pi 

-8^1  _  ^21 
2  2  * 

5p2  p2 


(2.11) 


(2.12) 


£[/?12]  —  ^^{^12  ^ 
n=\ 

c<  7;^)(  7^^  ><■  - 

n=l 

ji  -e)p^iP2  +  c) 

0p,p2  +  c[pi  +p2]+C^  ’ 

oo 

n=l 

oo 

(l-g)PlP2 

0P)P2  +  C[pj  +  P2]  +  ’ 


OO 

Pl%2  +  <^) 
0PlP2  +  c[p,  +  P2]  + 


(2.13) 


(2.14) 


(2.15) 
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where  c*'  is  the  mean  of  the  exponential  censoring  time.  Thus  for  N  individuals  the 
asymptotic  variances  of  the  estimators  are 


varm  =  m~'-, 

VarCpJ  =  Iip,r\ 

=  /(P2)~‘ 

w'here 

/?,o  Rn 

m  =  N  £[-f  + - 

6^  {\-df 

R\(\  4"  R\2 

/(p.)  =  iV  £[  -'°  ,  ]; 

Pi 

/(P2)  =  A'  ]. 

Pi 


(2.16) 

(2.17) 

(2.18) 

(2.19) 


The  expression  for  the  survival  function  S{t)  =  P{D  >  f)  for  this  continuous  time 
Markov  chain  is 


Sir)  =  { 


/i  +  Pi 

- ; -  g 

>•2 


Pi 

>1 


e^'']e 


Pi 

■^-1  “  ^-1 


(2.20) 


where  ?.j  are  the  roots  of  the  equation 

^PiP2+y(Pi +P2)+J'^  =  0-  (2-21) 

A 

The  maximum  likelihood  estimator,  denoted  as  for  the  survival 

probability  is  [Ref.  5:  p.  5  eqn  1.17] 


Pm{D>i}  = 


(2.22) 


where  A,  and  /j  are  roots  of  the  equation 
^p\Pi  +  ^(Pi  +  P2)  +  =  0. 
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(2.23) 


1 


Since  the  maximum  likelihood  estimators  are  orthogonal,  the  asymptotic  variance  of 
Pm{D  >  /}  [Ref.  5:  p.  5  eqn  1.19]  is  approximately 


VariPM{D> 


=  Varies 

oy  dpi  dp2 


where 


dS_S  dS  ^>-1  dS  d).2 
dd  d  5;.,  dd  cQ  ' 


es  S  dS  o>-\  8S  ^>1 


cpi  Pi  oAi  8p^  8A2  op^  ’ 


dp2 

—  A2)  >-2 

.  e'2‘  _ 

dS 

...-5  , 

0p, 

(3/, 

(;.i-;.2)  (/ 

1  ~  ‘^■2) 

cS 

5  , 

0p, 

a/2 

(>..->■2)  (>• 

1  ~  '^•2) 

dAi 

-  P1P2 

86 

Pi  +  P2  +  2/i, 

8;., 

—  Ai  —  6p2 

dpj 

Pi  +  P2  +  2/< 

8;., 

-  -  0Pi 

‘•1  ■'  ^P2 


-  ^'j  h 


^ ,  (fL. 

'‘•2 


OPI  Pi  +  P2  +  2;.,  ■ 


(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 


3.  Asymptotic  Renewal  Estimator 

In  this  subsection,  we  describe  an  asymptotic  renewal  estimator(A.R.E.)  for 
P{D  >  /}.  A  conditioning  argument  yields  the  following  equation  for  P{D  >  r}  ; 
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P{D  >  r)  =  g{i)  +  (1  -  0)  f  (ds)  P[D>t-  5}  (2.33) 

*^o 

with 

g{t)  =  F,{t)  +  (1  -  f  V,  ids)  F,it  -  s).  (2.34) 

*^o 

Thus,  P[D>  t)  satisfies  a  renewal- type  equation  with  defective  inter- renewal  distribution 

L(r)  =  (l-0)(f,*f2)(r)  (2.35) 

where  F^  is  the  sojourn  time  distribution  in  state  i,  F{t)  =  1  —  f;(r),  and  {F*F^[t)  denotes 
the  convolution  of  F^  and  fj.  Following  Feller  [Ref  6],  let  k  be  such  that 


Then,  under  certain  integrability  conditions,  if  {F*F^  is  not  arithmetic 
lim 

r-»oo  ^ 

where 


(2.37) 


(2.38) 


(2.39) 
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=  I  g{s)ds 
0 


oo 

e'“[F,(5)  +  (l-0)J  F,idu)F^is-u)ys 


= f  - 1} + (1  -  0)  umuk)  - 1}]. 


(2.40) 


Since 


it  follows  that 


+  ^  [«*)-!] 


1 

4>2ik) 


] 


=  X  [0iW-i]+ J  [1-(1-«)<^,W] 


=  y  [0. W  -  1]  +  y  [1  -  <t>M)  +  e  <t>,{k)-] 

=  y0.W.  (2.41) 

A  A 

Let  F;  be  the  Kaplan-Meier  estimate  of  F^,  and  6  be  the  maximum  likelihood 
estimate  of  6  [Ref.  5:  p.  10);  put 


J'*0O 

e^^F^ids). 

0 


(2.42) 
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The  asymptotic  renewal  estimator  (P^{Z)>r})  [Ref.  5:  p.  11  eqn  3.11]  of  the  survival 
probability  P{D>  t]  is 


PAD>t}  =  e^‘4- 


where  k  is  the  solution  to  the  equation 

(1-0)0,^#)  =  1; 


(2.43) 


(2.44) 


A 


(*oo 


(1-0) 


/^(f,*Fj)(rfs); 


(2.45) 


A  fi  f\  A 

b  =  (2.46) 

k 

A 

In  the  simulation  k  is  obtained  by  using  the  golden  section  search  method. 

The  asymptotic  renewal  estimator  is  undefined  if  all  the  sojourn  times  for  a 
particular  state  are  censored  since  the  Kaplan-Meier  estimator  is  not  defined  in  this  case. 

C.  CONFIDENCE  INTERVAL  PROCEDURES 

A  confidence  interval  for  an  unknowm  parameter  gives  both  an  indication  of  the 
numerical  value  of  the  unknown  parameter  and  a  measure  of  how  confident  we  are  of 
that  numerical  value.  Two  statistics  L  and  U  form  a  (1  —  a)  100%  confidence  interval  for 
0,  if  under  repeated  random  sampling  L^0<U  (1  —  a)  100%  of  the  time.  A  confi¬ 
dence  interval  procedure  for  an  unknown  parameter  0  is  also  used  to  make  a  decision 
concerning  0,  as  in  classical  hypothesis  testing  or  decision  making,  or  to  indicate  the 

A 

accuracy  and  variability  of  a  point  estimator  0. 

Confidence  interval  procedures  for  the  three  estimators  for  P{D>  t}  will  be  de¬ 
scribed  in  this  section,  starting  with  the  confidence  interval  procedure  for  the  Kaplan- 
Meier  estimator(K.M.E.).  Procedures  for  the  maximum  likelihood  estimator(M.L.E.) 
and  the  asymptotic  renewal  estimator  will  then  be  discussed. 

A  preliminary  transformation  to  approximately  symmetrize  the  sampling  distrib¬ 
ution  of  the  estimator  is  often  beneficial  [Ref  7).  For  this  study  we  consider  two  tran¬ 
sformations,  the  arc-sine  and  log  transformations,  w’hich  tend  to  stabilize  and  also 
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approximately  symmetrize  the  data.  These  transformations  were  suggested  by  the  work 
of  Gaver  and  Miller  [Ref  S).  Since  the  individuals  are  independent,  the  number  of  in¬ 
dividuals  surviving  a  fixed  time  t  would  have  a  binomial  distribution  if  there  were  no 
censoring.  The  logarithmic  and  arc-sine  transformations  have  been  beneficial  in  this 
case. 

The  confidence  intervals  for  the  Kaplan-Meier  estimator  and  the  maximum  likeli¬ 
hood  estimator  are  the  asymptotic  normal  confidence  intervals  using  the  transformed 
estimator. 

The  confidence  interval  for  the  arc-sine  transformed  estimator  is  computed  as  fol¬ 
lows.  Let  S(/)  be  the  estimator  of  S(/)  =  P{D  >  /}  .  Since 


d  ■  -1  r- 
— sm  jx 
dx 


1  1  1 

2  V-  ’ 


(2.47) 


a  Taylor  expansion  yields 

sin-'(x/-S(/) )  ^  sin-'(v'^)  +  -'t  ---  ..=^(5(0  -  5(r)).  (2.48) 

2  -  ■^W)S(0 

Hence,  the  approximate  variance  of  sin"'(>^/ S(r) )  is 

Far[sin-'(VS(/))]  ^  -j— — -  ?^«r[5(0].  (2.49) 

^  (l-5(r))5(/) 

IfS(f)  =  I,  we  set  Var[_  S{i) )]  equal  to  0.  Therefore,  a  (1  —  a)  100%  confidence 

interval  for  sin'\/S(r)  is 

sin-’(^/^)  ±  [z,_,/2]VLar[sin-'(VS(0)]  (2.50) 


where  is  the  (1  —  a/2)100%  point  of  a  standard  normal. 

A  confidence  interval  for  the  log  transformed  estimator  of  S(t)  is  computed  as  fol¬ 
lows.  Since 


-£-ln  X  =  4--  (2-51) 

a  Taylor  expansion  yields 

In^t)  ^  lnS{t)  +  ^(5(i)-S(0).  (2.52) 
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Thus  the  approximate  variance  of  In  S(r)  is 

Far[/«5(/)]  ^  (2.53) 

S(0' 

If  S(r)  =  1,  we  set  Far[S(/)]  =  0  .  A  (1  —  a)100%  confidence  interval  for  In  S{t)  is 

/«S(r)  ±  (2.54) 

where  is  the  (1  —  a/2)100%  point  of  a  standard  normal. 

1.  Asymptotic  Normal  Confidence  Intervals  for  K.M.E. 

A 

The  asymptotic  variance  of  the  Kaplan-Meier  estimator  (Pjf{D>/})  is  given 
by  equation  (2.2).  The  asymptotic  normal  confidence  intervals  for 
sin-'yj Pg{D  >  ;}  and  In  Pjf{D  >  /}  are  evaluated  using  equations  (2.50)  and  (2.54)  re¬ 
spectively.  The  corresponding  confidence  intervals  for  Pjf{Z)  >  r)  are  formed  by  inverse 
transformation.  If  the  lower  limit  is  less  than  0  it  is  set  equal  to  0  and  if  the  upper  limit 
is  greater  than  1  it  is  set  equal  to  1. 

2.  Asymptotic  Normal  Confidence  Intervals  for  M.L.E. 

The  asymptotic  variance  for  the  maximum  likelihood  estimator  (P*f(Z>  >  t})  is 
given  by  equation  (2.22).  The  asymptotic  normal  confidence  intervals  for 
sin-'(V  >  t} )  and  In  P*f{Z)  >  r}  can  be  constructed  using  equations  (2.50)  and 
(2.54)  respectively.  Confidence  intervals  for  P^{D>r}  can  be  obtained  by  inverse 
transformation.  If  the  resulting  confidence  interval  has  a  lower  limit  less  than  0  it  is  set 
equal  to  0  and  if  it  has  an  upper  limit  greater  than  1  it  is  set  equal  to  1. 

3.  Jackknife  Procedure  for  A.R.E. 

The  jackknife  technique  was  first  introduced  by  Quenouille  [Ref  9]  and  later 
utilized  by  Tukey  (Ref.  10]  for  bias  reduction  and  robust  interval  estimation.  A  review 
can  be  found  in  Miller  [Ref.  1 1).  The  jackknife  is  designed  to  do  various  jobs  fairly  well, 
however  it  is  desirable  to  avoid  (in  jackknifing)  sampling  distributions  with  (i)  abrupt 
ends  and  (ii)  one  or  more  straggling  tails,  and  it  is  probably  desirable  to  avoid  those  that 
are  strongly  unsymmetrical  [Ref  12].  Confidence  intervals  for  the  asymptotic  renewal 
estimator  will  be  obtained  using  the  jackknife  procedure  on  the  arc-sine  and  log  trans¬ 
formed  estimates. 

The  jackknife  procedure  for  the  confidence  interval  of  the  asymptotic  renewal 
estimator  (P^{D  >  r})  is  implemented  as  follows. 
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1.  Generate  data  for  N  individuals. 


A 

2.  Compute  >  /}  using  all  data. 

3.  Transform  P^{D  >  /}  into  In  P^{D  >  /}  and  sm~' yj P^{D  >  r}  which  will  be  denoted 

as  TL„„  and  TS,„  respectively.  Divide  the  N  individuals  into  n  subgroups  such 

N 

that  each  subgroup  contains  —  individuals. 

4.  Compute  >  r}  leaving  out  all  data  of  the  /'*  subgroup  and  transform  it  into 

A  /a 

In  Pj,{D>  t)  and  sin-',y  P^{D>  i)  which  will  be  denoted  as^,,  and  ys^  respectively. 

5.  Compute  the  pseudo-values(denoted  as  TL.,,  and  TS., ); 


Y L*]  —  n  Y  Ljji 

-  {n-  \)yli. 

(2.55) 

YS,i  =  nYS^i 

-  (n-  1) 

(2.56) 

6.  Compute  the  average  of  pseudo-values  which  is  the  jackknifed  estimate  for  the 
transformed  asymptotic  renewal  estimator(denoted  as  YL,  and  YS. ); 


-h 

# 

1 _ 1 

II 

...  +  yl.„i 

(2.57) 

YS,  =  + 

...  +  YS,„l 

(2.58) 

and 

compute  the  variance 

of  the  average 

of  the  pseudo-values(denoted  as 

vu, 

and  ra); 
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n 

7.  Compute  the  (approximate)  two-sided  (1  —  a)100%  confidence  intervals  for  each 
of  the  transformed  estimators  as  follows; 

YL,  ±  (2.61) 

and 

yS,  ±  (2-62) 

where  —  1)  is  the  (1  —  a/2)l00%  point  of  student's  t  with  n-1  degree  of  free¬ 
dom. 

8.  Confidence  intervals  for  P{D>  /}  are  obtained  by  inverse  transformation.  If  the 
resulting  interval  has  a  lower  limit  less  than  0  it  is  set  equal  to  0  and  if  it  has  an 
upper  limit  greater  than  1  it  is  set  equal  to  1. 
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III.  ANALYSIS  OF  SIMULATION  RESULTS  FOR  CONFIDENCE 

INTERVAL  PROCEDURES 

A.  SIMULATION 

All  simulations  were  carried  out  on  an  IBM  3033AP  computer  at  the  Naval  Post¬ 
graduate  School  using  the  LLRANDOM  II  random  number  generator  package[Ref. 
13].  Plots  of  simulated  estimates  and  confidence  intervals  were  produced  by  an  exper¬ 
imental  APL  package  GRAFSTAT  which  the  Naval  Postgraduate  School  is  using  under 
a  test  agreement  with  IBM  Watson  Research  Center,  Yorktown,  Heights,  NY. 

The  data  for  the  simulation  experiments  are  generated  as  follows; 

1 .  An  individual  starts  in  state  1  at  time  0. 

2.  An  exponential  censoring  time  with  mean  1/c  is  generated. 

3.  An  exponential  sojourn  time  in  state  1  with  mean  1/  p,  is  generated. 

4.  The  sojourn  time  and  censoring  time  are  compared;  if  the  sojourn  time  is  smaller, 
then  the  sojourn  time  is  recorded  and  given  an  uncensored  index  'O';  if  the  censor¬ 
ing  time  is  smaller,  the  sojourn  time  truncated  at  the  censoring  time  is  recorded  and 
given  a  censored  index  '1'.  In  the  latter  case,  the  death  time  is  recorded  as  the 
truncated  sojourn  time  and  associated  with  a  censored  index  of 'T  and  the  simu¬ 
lation  for  the  first  individual  is  completed. 

5.  If  the  process  is  not  censored  in  state  I,  a  uniform  random  number  is  generated 
and  compared  with  6;  if  less  than  0  the  process  jumps  to  state  0,  and  the  uncen¬ 
sored  death  time  with  index  '0'  is  recorded  and  the  simulation  for  the  first  individual 
is  completed;  if  greater  than  6,  the  process  jumps  to  state  2. 
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6.  An  exponential  sojourn  time  for  state  2  uith  mean  1/  Pj  is  generated  and  the  total 
time( sojourn  time  in  state  1  plus  sojourn  time  in  state  2)  is  compared  to  the  cen¬ 
soring  time  with  the  same  actions  as  listed  above. 

7.  If  the  process  is  not  censored  by  the  end  of  the  sojourn  time  in  state  2,  the  process 
jumps  to  state  1  and  continues  until  an  uncensored  or  censored  death  occurs.  The 
time  is  recorded  and  the  next  individual  is  started. 

This  procedure  continues  until  N  individuals'  data  have  been  generated.  Using  this 

A  A 

data,  the  Kaplan-Meier  estimate  (P^r))  ,  maximum  likelihood  estimate  (PM), 

A 

asymptotic  renewal  estimate  and  their  respective  90®/'o  and  80%  two-sided  confi¬ 

dence  intervals  are  computed.  This  completes  the  one  super-replication.The  simulation 
is  replicated  for  SR  =  500  super-replications  utilizing  dilTerent  seeds  to  generate  the 
data.  For  the  simulated  model  described  above  parameter  values  of 
Pi  =  U  Pj  =  U  ®  =  0.5,  and  c  =  ,  and  are  used,  and  for  each  super¬ 

replication  \  =  50  individuals'  data  are  generated. 

B.  ANALYSIS 

In  this  section  results  from  the  simulation  experiments  will  be  reported.  In  .Ap¬ 
pendix  A,  the  true  survival  probability,  w'hich  is  obtained  using  equation  (2.20),  is  given 
at  the  various  values  of  t.  The  times  considered  are  t  =  1.0,  3.0,  5.0,  7.0,  9.0,  11.0,  13.0, 
and  15.0. 

Tables  2  through  5  of  Appendix  B  report  the  number  of  super-replications  that  still 
have  defmed  Kaplan-.Meier  estimates  at  various  times  t.  The  column  headed  by  KMO 
reports  the  number  of  super-replications  that  have  defined  Kaplan-Meier  estimates  at 
time  t  for  the  death  time  survival  function.  The  columns  headed  KMl,  KM2,  repectively 
show  the  number  of  super-replications  still  defined  for  the  Kaplan-Meier  estimate  of  the 
survival  function  at  time  t  for  the  sojourn  times  in  state  1,  respectively  state  2.  These 
numbers  indicate  the  effect  the  censoring  has  on  the  Kaplan-Meier  estimate.  .As  ex¬ 
pected,  increasing  the  mean  time  to  censoring  increases  the  number  of  super-replications 
censored.  In  all  cases  the  survival  functions  for  the  sojourn  times  in  state  1  and  2  are  less 
heavily  censored  than  that  of  the  death  times. 

Plots  for  comparing  the  methods  investigated  for  the  Kaplan-Meier  estimator  are 
given  in  Figures  1  through  3.  The  difference  between  the  mean  of  the  estimated  survival 

A 

probability  and  the  true  survival  probability  (S(f)  —  P[D  >  /}),  the  relative  differences 
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((S(/)  —  P{Z)  >/})/?{£)> /}),  and  the  relative  root  mean  square  errors 
(RMSE/P{D  >  /})  are  respectively  plotted  in  Figures  1  through  3.  The  numerical  results 
are  recorded  in  Tables  6  through  17.  The  means,  and  root  mean  square  errors{RMSE) 
from  SR  =  500  super- replications  are  computed  as  follows; 


SR 


1=1 


and 


RMSE 


SR 


'-TrT. 


1=1 


(3.2) 


where  .S,(/)  is  the  point  estimate  of  the  true  value  S(/)  at  time  t  in  the  i’’'  super-replication 
and  SR  is  the  number  of  super-replication. 

If  a  Kaplan-Meier  estimate  is  undefined,  its  value  is  taken  to  be  its  last  defined 
value.  The  modifications  to  the  Kaplan-Meier  estimate,  MOD  1  and  MOD  2,  are  de¬ 
scribed  in  Chapter  2.  Since  the  unmodified  Kaplan-Meier  survival  function  may  not 
tend  to  0  as  t  -»  oo  ,  it  tends  to  overestimate  the  true  survival  probability.  On  the  other 
hand,  since  the  modifications  MOD  1  and  MOD  2  make  the  undefined  Kaplan-Meier 
sur%'ival  function  honest,  they  do  well  for  large  t  and  have  less  variability  than  the  un¬ 
modified  K.M.E..  However  the  modifications  appear  to  bias  the  estimates  for  moderate 
times  t.  All  these  methods  improve  as  the  mean  censoring  time  decreases.  MOD  1  does 
slightly  better  than  MOD  2. 

In  Figures  4  through  9,  plots  are  presented  for  the  comparison  between  the  three 
estimators.  The  differences,  relative  differences,  and  relative  RMSE's  are  computed  and 
plotted  in  the  same  manner  as  before.  Figures  4  through  6  show  results  for  the 
Kaplan-Meier  and  asymptotic  renewal  estimators  using  the  u  nmodified  Kaplan-Meier 
estimator.  Figures  7  through  9  show  the  results  obtained  when  the  Kaplan-Meier  esti¬ 
mates  are  modified  to  be  MOD  1  in  both  the  K.M.E.  and  A.R.E..  The  numerical  results 
are  recorded  in  Appendix  C.  As  expected,  the  maximum  likelihood  estimator(M.L.E.), 
which  uses  the  most  correct  information  about  the  model,  tends  to  have  the  smallest 
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relative  RMSE's  and  difTerences  between  estimated  and  true  means.  In  the  case  of 
greatest  censoring,  the  M.L.E.  shows  a  slight  bias. 

The  Kaplan-Meier  estimator(K.M.E.),  which  uses  the  least  information  about  the 
model,  does  well  for  small  times  even  in  the  greatest  censoring  case,  but  it  is  worse  than 
the  others  for  large  times  due  to  the  undefined  estimates.  It  has  the  largest  values  of 
mean  square  error  among  the  three  estimators.  Theiefore  the  K.M.E.  tends  to  be  more 
variable  than  the  other  estimators. 

Not  surprisingly,  the  asymptotic  renewal  estimator(A.R.E.),  which  uses  an 
asymptotic  exponential  approximation  to  the  probability,  has  high  relative  RMSE's  and 
differences  between  estimated  and  true  means  for  small  times  t.  As  t  increases  the  rela¬ 
tive  RMSE's  and  differences  between  the  means  decrease.  Tables  6  through  9  indicate 
that  for  moderate  to  large  times  the  A.R.E.  has  RMSE's  that  are  comparable  or  less 
than  those  of  the  M.L.E..  The  time  at  which  they  become  comparable  is  a  function  of 
the  amount  of  censoring.  If  the  mean  censoring  time  is  larger,  the  R.MSE's  become 
comparable  sooner.  Tables  6  through  12  indicate  that  while  modifying  the  Kaplan-Meier 
estimate  of  the  survival  function  of  the  death  times  improves  it  for  large  times  t,  it  cre¬ 
ates  a  bias  for  moderate  times  t.  Using  modified  Kaplan-Meier  estimators  of  the  surv’ival 
functions  for  the  sojourn  times  in  states  1  and  2  does  not  improve  the  asymptotic  re¬ 
newal  estimator. 

In  Figures  10  through  17,  plots  are  presented  for  comparing  the  confidence  interval 
procedures  of  the  three  estimators.  In  order  to  compare  the  performance  of  the  confi¬ 
dence  interval  procedures,  we  use  the  following  measures:  the  coverage  fraction(CVR), 
the  average  half  length(AHL)  of  the  confidence  interval,  and  the  standard 
deviation(SHL)  of  the  half  lengths.  The  AHL  is  determined  by  summing  the  half  lengths 
of  the  confidence  intervals  of  the  replications  and  dividing  that  sum  by  the  number  of 
replications(SR).  The  SHE  is  computed  by  summing  the  square  of  the  differences  be¬ 
tween  the  individual  half  lengths  and  the  average  half  length,  dividing  that  value  by 
(SR-1)  and  taking  the  square  root  of  the  resulting  value.  Among  confidence  interval 
procedures  which  achieve  the  desired  coverage  rate(0.90,  and  0.80),  the  confidence  in¬ 
terval  procedure  which  yields  the  smallest  AHL  is  preferred.  Also  preferred  is  a  small 
SHL  representing  a  stable  confidence  interval  procedure. 

For  each  procedure,  the  number  of  intervals  covering  the  true  value  P{D  >  r}  is 
recorded  as  well  as  the  number  of  intervals  that  are  too  high  or  too  low.  These  results 
are  reported  in  Tables  18,  20,  22,  24,  26,  28,  30,  and  32  (Appendix  D).  Next  to  each 
coverage  count  is  given  the  corresponding  coverage  fraction  in  parentheses.  If  a 
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(I  —  a)  100%  confidence  interval  procedure  is  performing  well,  then  this  interval  should 
cover  about  (1  —  a)  100%  of  time.  A  95%  confidence  interval  for  the  coverage  fraction 
is  computed  using  P  ±  1.96lP(l  —  P)/5/?P  .where  P  is  the  proportion  of  (1  -  a)%  con¬ 
fidence  intervals  that  cover  the  true  value  of  P(D  >  /}  and  SR  is  the  number  of  super¬ 
replications.  Therefore,  if  a  80%(90%)  confidence  interval  procedure  is  working  well, 
then  out  of  500  super-replications  between  382(436)  and  418(463)  confidence  intervals 
should  cover  the  true  value.  These  intervals  correspond  to  0.7649(0.8351)  and 
0.8737(0.9263)  respectively.  The  average  half  lengths  of  the  confidence  intervals  for 
S(r)  are  reported  in  tables  19,  21,  23,  25,  27,  29,  31,  and  33  (Appendix  D).  The  standard 
deviation  of  the  half  length  is  given  in  parentheses  below  the  average  half  length.  If  an 
estimator  is  performing  well,  its  confidence  interval  should  not  only  have  the  correct 
coverage  fraction  but  also  a  small  average  half  length. 

In  Figures  10  through  13,  the  coverage  fraction  and  the  average  half  length  of 
confidence  intervals  using  the  log  transformation  for  each  estimator  are  presented.  In 
Figures  14  through  17,  the  coverage  fraction  and  the  average  half  length  of  confidence 
intervals  using  the  arc-sine  transformation  are  presented.  The  two  horizontal  lines  in  the 
coverage  plot  show  the  95%  confidence  interval  for  the  coverage  fraction. 

The  following  remarks  concern  the  confidence  intervals  obtained  using  the  arc-sine 
transformation.  The  asymptotic  normal  confidence  intervals  for  the  M.L.E.  have  the 
correct  coverage  for  small  to  moderate  times  but  tend  to  slightly  undercover  for  large 
times.  The  asymptotic  nomal  confidence  intervals  for  the  K.M.E.  have  the  correct  cov¬ 
erage  for  small  times  but  undercover  for  moderate  to  large  times  due  to  the  undefined 
estimates.  As  the  mean  censoring  time  increases,  the  coverage  of  the  Kaplan-Meier 
confidence  intervals  improves.  The  confidence  intervals  for  the  asymptotic  renewal  esti¬ 
mate  tend  to  undercover.  Using  a  modified  Kaplan-Meier  estimate  makes  very  little  dif¬ 
ference  in  the  results. 

The  following  remarks  concern  the  confidence  intervals  obtained  by  using  the  log 
transformation.  Once  again  the  confidence  intervals  for  the  K.M.E.  undercover  for 
moderate  and  large  times  t.  The  coverage  is  slightly  worse  than  that  obtained  using  the 
arc-sine  transformation.  The  confidence  intervals  for  the  M.L.E.  have  the  correct  cov¬ 
erage  and  the  smallest  AHL.  The  confidence  intervals  for  the  Jackknifed  A.R.E.  have 
the  correct  coverage  for  all  but  the  80%  interval  at  t  =  1  for  c=  1/2.  The  average  half 
lengths  of  the  A.R.E.  confidence  intervals  are  larger  than  those  of  the  M.L.E.  confidence 
intervals.  Using  a  modified  Kaplan-Meier  estimate  makes  very  little  difierence  in  the  re¬ 
sults. 
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DIFFERENCE  BETWEEN  MEAN  AMD  TRUE  VALUE 


TIME  TIME 


RELATIVE  RMSE 


^nyi-win 


3nai-NV3w 


COMPARISON  BETWEEN  ESTIMATORS(UNMOD): 
between  Mean  and  True  value. 
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COMPARISON  BETWEEN  ESTIMATORS(UNMOD):  Relative 
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between  Mean  and  True  value. 
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RELATIVE  DIFFERENCE(MOD  1) 
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(C  -  1/2)  UNMOD  COVERAGE  (C  -  1/10) 


01  BO  90  to  ZO  0  Zl  B  ^  0 

3Svy3A03  anal-;-  ihv 


Figure  12.  CONFIDENCE  INTERVAL(90%,  LOG,  MOD  I):  Coverage 
Fraction  and  Average  Half  Length. 
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80  PERCENT  C.l  (LOG) 

(C  -  1/2)  MOD  1  COVERAGE  (C  -  1/10) 


Figure  13.  CONFIDENCE  INTERVAL(80%,  LOG,  MOD  1):  Coverage 
Fraction  and  Average  Half  Length. 
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PERCENT  C.l  (ARC-SINE) 

UNMOD  COVERAGE  (C  -  1/10) 


PERCENT  C.l  (ARC-SINE) 

UN  MOD  COVERAGE  (C  -  1/10) 


90  PERCENT  C.l  (ARC-SINE) 

(c  -  1/2)  MOD  1  COVERAGE  (C  .  1/10) 
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IV.  CONCLUSIONS 


This  thesis  considers  the  problem  of  estimating  the  survival  probability  P{D>  t)  for 
the  first  passage  time  to  state  0  for  a  semi-Markov  process  using  censored  data.  Simu¬ 
lation  is  used  to  study  the  small  sample  behavior  of  three  estimators  and  their  confidence 
interval  procedures. 

One  of  the  estimators  studied  is  the  Kaplan-Meier  estimator  of  the  first  passage 
times  to  state  0.  Both  the  unmodified  Kaplan-Meier  estimator  and  two  modifications, 
MOD  1  and  .MOD  2  making  the  estimated  distribution  honest  are  considered.  Another 
estimator  is  the  maximum  likelihood  estimator.  A  third  estimator,  the  asymptotic  re¬ 
newal  estimator,  uses  an  exponential  approximation  to  the  survival  function. 

The  following  conclusions  are  drawn  from  the  simulation  experiment. 

1.  The  modified  Kaplan-Meier  estimators  MOD  1  and  MOD  2  using  the  first  passage 
times  to  state  0  have  a  smaller  bias  for  large  times  than  the  unmodified  K.M.E.. 
However,  in  the  medium  range  of  times,  the  two  modified  procedures  MOD  1  and 
MOD  2  have  a  larger  bias  than  that  of  the  unmodified  K.M.E..  MOD  1  is  slightly 
better  than  MOD  2. 

2.  Modifying  the  Kaplan-Meier  estimates  of  the  sojourn  time  distributions  in  the 
asymptotic  renewal  estimate  does  not  improve  its  performance. 

3.  The  asymptotic  normal  confidence  intervals  for  the  Kaplan-Meier  estimator  of  the 
first  passage  times  to  state  0  using  the  arc-sine  transformation  have  a  slightly  better 
coverage  than  those  using  the  log  transformation. 

4.  The  confidence  intervals  using  the  log  transformed  estimators  are  preferred  to 
those  using  the  arc-sine  transformed  estimators  for  the  maximum  likelihood  esti¬ 
mator  and  the  asymptotic  renewal  estimator. 
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5.  The  confidence  intervals  for  the  jackknifed  log  transformed  asymptotic  renewal 
estimator  have  the  correct  coverage  for  all  but  the  smallest  times.  This  estimator 
makes  no  assumptions  concerning  the  parametric  form  of  the  sojourn  time  distrib¬ 
utions.  It's  expected  that  it  will  also  perform  well  in  cases  in  which  the  sojourn  time 
distributions  are  not  exponential. 

6.  The  confidence  intervals  for  the  maximum  likelihood  estimator  have  the  correct 
coverage,  and  also  have  the  smallest  average  half  length.  However  the  estimators 
depend  on  the  parametric  form  of  the  sojourn  time  distributions.  If  the  parametric 
form  is  incorrectly  specified,  the  M.L.E.  can  be  quite  biased,  Kim[Rer  1]. 
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APPENDIX  A.  TRUE  SURVIVAL  PROBABILITY 


Table  1.  TRUE  SURVIVAL  PROBABIFITY  FOR  MODEL 


TIME 

P{D  >  /} 

1.00 

0.6644 

3.00 

0.3554 

5.00 

0.1974 

7.00 

0.1099 

9.00 

0.0612 

11.00 

0.0340 

13.00 

0.0189 

15.00 

0.0105 
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APPENDIX  B.  NUMBER  OF  K.M.E.  DEFINED  AT  TIME  T 


Table  2.  NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 
C=  1/2 


TIME 

KMO 

KMl 

KM2 

1.00 

500 

500 

496 

496 

428 

369 

366 

314 

334 

7.00 

215 

308 

333 

9.00 

171 

308 

333 

11.00 

165 

308 

333 

13.00 

164 

308 

333 

15.00 

164 

308 

333 

Table  3.  NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 
C=  1/5 


TIME 

KMO 

KMl 

KM2 

1.00 

500 

500 

500 

500 

488 

463 

496 

436 

425 

7.00 

456 

421 

419 

9.00 

375 

420 

419 

11,00 

330 

420 

419 

13.00 

309 

420 

419 

15.00 

305 

420 

419 

Table  4.  NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 
C=  1/10 


TIME 

KMO 

KMl 

KM2 

1.00 

500 

500 

497 

484 

499 

462 

455 

7.00 

489 

451 

451 

9.00 

466 

447 

450 

11.00 

436 

447 

450 

13.00 

412 

447 

450 

15.00 

399 

447 

450 

Table  5.  NUMBER  OF  KAPLAN-MEIER  ESTIMATES  DEFINED  AT  TIME  T: 
C=  1/100 


TIME 

KMO 

KMl 

KM2 

1.00 

500 

500 

500 

499 

500 

498 

494 

7.00 

500 

494 

493 

9.00 

499 

494 

493 

11.00 

499 

494 

493 

13.00 

498 

494 

493 

15.00 

494 

494 

493 
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APPENDIX  C.  STATISTICS  OF  THREE  ESTIMATORS 


Table  6.  STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/2 


Table  7.  STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/5 
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Table  8.  STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/10 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6628 

0.0046 

0.6604 

0.6428 

0.0031 

mEm 

0.3518 

0.3535 

0.3476 

0.1962 

0.1981 

0.0026 

0.1905 

KiKh 

7.0 

0.1087 

0.1124 

0.0015 

0.1057 

9.0 

0.0644 

0.0643 

0.0008 

0.0593 

0.0008 

11.0 

0.0017 

0.0370 

0.0004 

0.0336 

0.0004 

13.0 

0.0249 

0.0013 

0.0215 

0.0002 

0.0192 

0.0002 

15.0 

0.0179 

0.0011 

0.0126 

0.0001 

0.0111 

0.0001 

Table  9.  STATISTICS  OF  THREE  ESTIMATORS  :  UNMOD,  C=  1/100 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6633 

0.0044 

0.6591 

3.0 

0.3558 

0.3513 

5.0 

0.1974 

0.1956 

WSmSM 

0.0021 

7.0 

0.1078 

0.1100 

0.0012 

0.0012 

9.0 

0.0593  - 

wnSuH 

0.0623 

0.0006 

0.0006 

11.0 

0.0327 

0.0007 

0.0355 

0.0003 

0.0347 

0.0003 

13.0 

0.0179 

0.0004 

0.0204 

0.0001 

0.0198 

0.0001 

15.0 

0.0110 

0.0003 

0.0118 

0.0001 

0.0114 

0.0001 
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Table  10.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/2 


Table  11.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/5 
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Table  12.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/10 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6628 

0.6604 

0.6247 

mSm 

0.3518 

0.3535 

0.3389 

0.0056 

HI 

0.1962 

0.1981 

0.1862 

0.0034 

7.0 

0.1074 

0.1124 

0.0015 

0.1035 

0.0018 

9.0 

0.0605 

0.0643 

0.0008 

0.0582 

0.0008 

11.0 

0.0304 

0.0016 

0.0370 

0.0004 

0.0330 

0.0004 

13.0 

0.0149 

0.0009 

0.0215 

0.0002 

0.0189 

0.0002 

15.0 

0.0060 

0.0004 

0.0126 

0.0001 

0.0109 

0.0001 

Table  13.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  1,  C=  1/100 


TIME 

K.M.E 

M.L.E 

A.R.E 

MEAN 

SRMSE 

MEAN 

SRMSE 

MEAN 

SRMSE 

1.0 

0.6633 

0.0044 

0.6591 

■m 

HI 

0.3558 

0.0045 

0.3513 

HI 

0.1974 

0.0034 

0.1956 

0.0021 

7.0 

0.1078 

0.0019 

0.1100 

0.1088 

0.0012 

9.0 

0.0592 

0.0011 

0.0623 

0.0612 

0.0006 

11.0 

0.0326 

0.0007 

0.0355 

0.0347 

0.0003 

13.0 

0.0178 

0.0004 

0.0204 

0.0001 

0.0198 

0.0001 

15.0 

0.0108 

0.0003 

0.0118 

0.0114 

0.0001 
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Table  14.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/2 


Table  15.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/5 
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Table  16.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C=  1/10 


Table  17.  STATISTICS  OF  THREE  ESTIMATORS  :  MOD  2,  C  =  1/ 100 
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APPENDIX  D.  CONFIDENCE  INTERVALS 


Table  18.  TWO-SIDED  90  %  COVERAGE  FRACTION(UNMOD,  C=  1/2) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

toohigh 

44(.09) 

30{.06) 

42(.08) 

26(.05) 

12(.02) 

11(.02) 

1.0 

cover 

444(.89) 

453(.91) 

445(.89) 

455(.91) 

434(.87) 

401(.80) 

too  low 

12(.02) 

17(.03) 

13(.03) 

19(.04) 

54(.ll) 

88(.18) 

too  high 

55(.ll) 

39{.08) 

35(.07) 

29(.06) 

12(.02) 

10(.02) 

3.0 

cover 

439(.8S) 

438{.88) 

449(.90) 

447(.89) 

453(.91) 

402(.80) 

too  low 

6(.01) 

23(.05) 

16(.03) 

24(.05) 

35{.07) 

88(.18) 

too  high 

100(.20) 

78(.16) 

39(.08) 

27{.05) 

12(.02) 

9(.02) 

5.0 

cover 

32I(.64) 

340(.68) 

444(.89) 

448(.90) 

452(.90) 

382(.76) 

too  low 

79(.16) 

82(.16) 

17(.03) 

25(.05) 

36(.07) 

109(.22) 

too  high 

162(.32) 

134(.27) 

41(.08) 

26(.05) 

12(.02) 

8(.02) 

7.0 

cover 

194{.39) 

222{.44) 

443(.89) 

448(.90) 

453(.91) 

355(.71) 

too  low 

144(.29) 

144(.29)  ; 

16(.03) 

26(.05) 

35(.07) 

137(.27) 

too  high 

224(.45) 

190(.38)  - 

43(.09) 

25(.05) 

12(.02) 

8(.02) 

9.0 

cover 

115(.23) 

149(.30) 

440(.88) 

440(.88) 

452(.90) 

344(.69) 

too  low 

161(.32) 

161C.32) 

I7(.03) 

35(.07) 

36{.07) 

148(.30) 

too  high 

285(.57) 

217(.43) 

46(.09) 

24{.05) 

12(.02) 

8(.02) 

11.0 

cover 

51(.10) 

119(.24) 

437{.87) 

435(.87) 

454{.91) 

322(.64) 

too  low 

164(.33) 

164(.33) 

17(.03) 

41(.08) 

34(.07) 

I70(.34) 

too  high 

323(.65) 

250(.50) 

46{.09) 

19(.04) 

12(.02) 

9{.02) 

13.0 

cover 

13(.03) 

86(.17) 

437(.87) 

437(.87) 

454(.91) 

305(.61) 

too  low 

164(.33) 

164(.33) 

17{.03) 

44(.09) 

34{.07) 

186(.37) 

too  high 

336(.67) 

290(.58) 

47(.09) 

16(.03) 

12(.02) 

9(.02) 

15.0 

cover 

46(.09) 

436(.87) 

434(.87) 

454(.91) 

285(.57) 

too  low 

164(.33) 

164(.33)  ! 

17(.03) 

50(.10) 

34(.07) 

206(.41) 
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Table  19.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  UNMOD,  C=  1/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.1245  .1223 

.0959  .0949 

.2032  .1862 

(.0107)  (.0102) 

(.0080)  (.0077) 

(.0841)  (.0719) 

3.0 

.1783  .1647 

.1337  .1289 

.2583  .1948 

(.0341)  (.0249) 

(.0108)  (.0108) 

(.1215)  (.0824) 

5.0 

.1775  .1494 

.1251  .1164 

.2572  .1490 

(.0904)  (.0722) 

(.0223)  (.0223) 

(.1506)  (.0840) 

7.0 

.1535  .1267 

.1045  .0927 

.2415  .1031 

(.1083)  (.0864) 

(.0295)  (.0287) 

(.1704)  (.0778) 

9.0 

.1463  .1207 

.0830  .0695 

.2243  .0683 

(.1108)  (.0887) 

(.0320)  (.0299) 

(.1831)  (.0660) 

11.0 

.1452  .1197 

.0644  .0505 

.2091  .0444 

(.1113)  (.0891) 

(.0316)  (.0281) 

(.1914)  (.0531) 

13.0 

.1452  .1197 

.0495  .0361 

.1966  .0285 

(.1113)  (.0891) 

(.0295)  (.0249) 

(.1969)  (.0414) 

15.0 

.1452  .1197 

.0379  .0256 

.1866  .0183 

(.1113)  (.0891) 

(.0267)  (.0213) 

(.2009)  (.0316) 
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Table  20.  TWO-SIDED  80  %  COVERAGE  FRACTION(UNMOD,  C=  1/2) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1 

too  high 

68(.14) 

59(.12) 

1  61(.12) 

54(.ll) 

33(.07) 

29(.06) 

1.0 

cover 

406(.81) 

396(.79) 

407(.81) 

405(.81) 

374(.75) 

342(.68) 

too  low 

26(.05) 

45(.09) 

32(.06) 

41(.08) 

93(.19) 

129(.26) 

too  high 

94(.19) 

75(.15) 

62(.12) 

47(.09) 

29(.06) 

26(.05) 

3.0 

cover 

381(.76) 

382(.76) 

404(.80) 

409(.82) 

401(.80) 

343(.69) 

too  low 

25(.05) 

43(.09) 

34(.07) 

44{.09) 

70(.I4) 

131(.26) 

too  high 

I29(.26) 

I08(.22) 

64(.13) 

48(.I0) 

25(.05) 

19{.04) 

5.0 

cover 

291(.58) 

302(.60) 

406{.81) 

405(.81) 

407(.81) 

328(.76) 

too  low 

80(.16) 

90(.18) 

30(.06) 

47(.09) 

68(.14) 

153(.31) 

too  high 

191(.38) 

164(.33) 

65(.13) 

48(.10) 

24(.05) 

17(.03) 

7.0 

cover 

165(.33) 

192{.38) 

406(.81) 

398(.80) 

403(.81) 

307(.61) 

too  low 

144(.29) 

144(.29) 

29(.06) 

54(.ll) 

73(.15) 

176(.35) 

too  high 

257(.51) 

217(.43) 

67(.13) 

47(.09) 

22(.04) 

16{.03) 

9.0 

cover 

82(.16) 

122{.24) 

405(.81) 

396(.79) 

403(.81) 

289(.58) 

too  low 

161(.32) 

161(.32) 

28(.06) 

57(.ll)  ■ 

75(.15) 

195(.39) 

too  high 

307(.61) 

267(.53)  ^ 

66(.13) 

45{.09)  ' 

21  (.04) 

15(.03) 

11.0 

cover 

29(.06) 

69(.14) 

405(.81) 

392(.78) 

403(.81) 

273(.55) 

too  low 

164(.33) 

164(.33) 

29(.06) 

63{.13) 

76{.15) 

212(.42) 

too  high 

334(.67) 

308(.62) 

69(.14) 

39(.08) 

21{.04) 

12(.02) 

13.0 

cover 

2(.00) 

28(.06) 

402(.80) 

396(.79) 

404(.81) 

267(.53) 

too  low 

I64(.33) 

I64(.33) 

29(.06) 

65(.13) 

75(.15) 

221  (.44) 

too  high 

336<.67) 

333(.67)  1 

70(.14) 

34<.07) 

21(.04)  1 

12(.02) 

15.0 

cover 

O(.OO) 

3(.01)  ; 

400(.80) 

400(.80) 

404(.81)  1 

252(.50) 

too  low 

164(.33) 

164(.33) 

30(.06) 

66(.13) 

75(.15) 

236(.47) 

Table  21.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  UNMOD,  C=  1/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.0968  .0957 

.0746  .0742 

.1596  .1447 

(.0082)  (.0080) 

(.0062)  (.0061) 

(.0748)  (.0602) 

3.0 

.1359  .1296 

.1032  .1010 

.2048  .1510 

(.0240)  (.0200) 

(.0085)  (.0085) 

(.1153)  (.0677) 

5.0 

.1311  .1181 

.0953  .0912 

.2010  .1152 

(.0655)  (.0574) 

(.0175)  (.0175) 

(.14300  (.0677) 

7.0 

.1126  .1002 

.0781  .0726 

.1845  .0795 

(.0784)  (.0686) 

(.0229)  (.0226) 

(.1598)  (.0621) 

9.0 

.1073  .0955 

.0607  .0545 

.1677  .0526 

(.0804)  (.0704) 

(.0245)  (.0236) 

(.1701)  (.0522) 

11.0 

.1065  .0947 

.0459  .0396 

.1534  .0340 

(.0808)  (.0707) 

(.0237)  (.0221) 

(.1765)  (.0417) 

13.0 

.1065  .0947 

.0343  .0283 

.1420  .0218 

(.0808)  (.0707) 

(.0217)  (.0196) 

(.1805)  (.0322) 

15.0 

.1065  .0947 

.0256  .0200 

.1330  .0140 

(.0808)  (.0707) 

(.0192)  (.0167) 

(.1832)  (.0245) 
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Table  22.  TWO-SIDED  90  %  COVERAGE  FRACTION(UNMOD,  C=  I/IO) 


xiMir 

COV¬ 

K.M.E 

M.L.E 

A.R.E 

i  IML 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

34(.07) 

26(.05) 

31(.06) 

24(.05) 

16(.03) 

12(.02) 

1.0 

cover 

453(.91) 

454(.91) 

448(.90) 

443(.89) 

436(.87) 

428(.86) 

too  low 

13(.03) 

20(.04) 

21(.04) 

33(.07) 

48(.10) 

60(.12) 

too  high 

24(.05) 

16(.03) 

27(.05) 

21(.04) 

24(.05) 

15(.03) 

3.0 

cover 

460(.92) 

456(.91) 

450(.90) 

444{.89) 

460(.92) 

450(.90) 

too  low 

16(.03) 

28(.06) 

23(.05) 

35(.07) 

16(.03) 

35(.07) 

too  high 

37(.07) 

20(.04) 

27(.05) 

19(.04) 

23(.05) 

13(.03) 

5.0 

cover 

453(.91) 

452(.90) 

452(.90) 

445(.89) 

457(.91) 

451(.90) 

too  low 

10(.02) 

28(.06) 

21(.04) 

36(.07) 

20(.04) 

36(.07) 

too  high 

37(.07) 

25(.05) 

27(.05) 

18(.04) 

23(.05) 

12(.02) 

7.0 

cover 

438(.88) 

434(.87) 

454(.91) 

442(.88) 

456(.91) 

444(.89) 

too  low 

25(.05) 

4J(.08) 

19(.04) 

40{.08) 

21{.04) 

44(.09) 

too  high 

47(.09) 

27(.05) 

26(.05) 

16<.03) 

24(.05) 

9(.02) 

9.0 

cover 

357(.71) 

377(.75)  ' 

455(.91) 

438(.88) 

455(.91) 

438(.88) 

too  low 

96(.19) 

96(.19)  ! 

19(.04) 

46(.09) 

21(.04) 

53(.ll) 

too  high 

59{.12) 

29(.06)  : 

26(.05) 

12(.02) 

24(.05) 

7(.01) 

11.0 

cover 

233(.47) 

263(.53) 

455{.91) 

434(.87) 

455(.91) 

428(.86) 

too  low 

208(.42) 

208(.42)  , 

19(.04) 

54(.ll) 

21(.04) 

65(.13) 

too  high 

59(.12) 

38(.08)  1 

26{.05) 

11(.02) 

24(.05) 

6(.01) 

13.0 

cover 

157(.31) 

175(.35) 

454(.91) 

433(.87) 

455(.91) 

422(.84) 

too  low 

287(.57) 

287(.57) 

20(.04) 

56(.ll) 

21  (.04) 

72(.14) 

too  high 

83(.17) 

35(.07)  1 

26(.05) 

7(.01) 

24(.05) 

4(.01) 

15.0 

cover 

73(.15) 

121(.24) 

456(.91) 

436(.87) 

455(.91) 

419(.84) 

too  low 

344{.69) 

344(.69) 

18{.04) 

57(.ll) 

21  (.04) 

77(.15) 

Table  23.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  UNMOD,  C=l/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.1121  .1105 

.0718  .0714 

.0932  .0903 

(.0064)  (.0062) 

(.0061)  (.0060) 

(.0317)  (.0279) 

3.0 

.1223  .1185 

.0972  .0953 

.1212  .1100 

(.0062)  (.0066) 

(.0042)  (.0044) 

(.0421)  (.0297) 

5.0 

.1123  .1046 

.0871  .0837 

.1160  .0969 

(.0132)  (.0142) 

(.0101)  (.0103) 

(.0522)  (.0299) 

7.0 

.0988  .0847 

.0688  .0644 

.0981  .0738 

(.0285)  (.0267) 

(.0134)  (.0135) 

(.0576)  (.0280) 

9.0 

.0806  .0640 

.0512  .0464 

.0789  .0523 

(.0437)  (.0363) 

(.0139)  (.0137) 

(.0596)  (.0241) 

11.0 

.0573  .0432 

.0369  .0321 

.0624  .0355 

(.0513)  (.0398) 

(.0128)  (.0122) 

(.0601)  (.0195) 

13.0 

.0416  .0305 

.0261  .0217 

.0494  .0235 

(.0506)  (.0380) 

(.0109)  (.0101) 

(.0601)  (.0151) 

15.0 

.0306  .0223 

.0182  .0145 

.0395  .0153 

(.0472)  (.0350) 

(.0089)  (.0079) 

(.0602)  (.0113) 
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Table  24.  TWO-SIDED  80  %  COVERAGE  FRACTION(UNMOD,  C=  I/IO) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

6I(.12) 

47(.09) 

^  56{.ll) 

48(.10) 

38(.08) 

31  (.06) 

1.0 

cover 

404(.81) 

402(.80) 

400(.80) 

403(.81) 

382(.76) 

379(.76) 

too  low 

35(.07) 

51(.10) 

44(.09) 

49(.10) 

80(.16) 

90(.18) 

too  high 

49(.10) 

38(.08) 

52(.10) 

42(.08) 

57(.ll) 

43(.09) 

3.0 

cover 

415(.83) 

417(.83) 

401(.80) 

403(.81) 

401(.80) 

398(.80) 

too  low 

36(.07) 

45(.09) 

47(.09) 

55(.ll) 

42(.08) 

59(.12) 

too  high 

56(.ll) 

47(.09) 

52(.10) 

40(.08) 

53(.ll) 

40(.08) 

5.0 

cover 

411(.82) 

392(.78) 

397(.79) 

395(.79) 

406(.81) 

395(.79) 

too  low 

33(.07) 

61(.12) 

51(.10) 

65(.13) 

41(.08) 

65(.13) 

too  high 

55(.ll) 

43{.09) 

51(.10) 

39(.08) 

50(.01) 

35(.07) 

7.0 

cover 

410(.82) 

390(.78)  1 

395(.79) 

396(.79) 

407(.81) 

392(.78) 

too  low 

35(.07) 

67(.13  ! 

54(.ll) 

65(.13) 

43(.09) 

73(.15) 

too  high 

7'(.14) 

53(.ll) 

51(.10) 

34(.07) 

50(.10) 

32(.06) 

9.0 

cover 

333(.67) 

351(.70) 

394(.79) 

399(.80) 

409(.82) 

387(.77) 

too  low 

96(.19) 

96(.19) 

55(.12) 

67(.13) 

41(.08) 

81(.16) 

too  high 

79(.16) 

59(.l2) 

53(.ll) 

31  (.06) 

49(.10) 

26(.05) 

11.0 

cover 

213(.43) 

233(.47) 

392(.78) 

400(.80) 

408(.82) 

3S2(.76) 

too  low 

208(.42) 

208(.42)  1 

55(11) 

69(.I4) 

43(.09) 

92(.18) 

too  high 

84(.17) 

56(.ll)  ■ 

55(.ll) 

29(.06) 

49(.10) 

25(.05) 

13.0 

cover 

129(.26) 

i57(.3i) ; 

391(.78) 

398(.89) 

406(.81) 

379(.76) 

too  low 

287(.57) 

287(.57) 

54(.ll) 

73(.15) 

45(.09) 

96(.19) 

too  high 

110(.22) 

58(.l2) 

54(.ll) 

26(.05) 

50(.I0) 

21(.04) 

15.0 

cover 

46(.09) 

98(.20) 

392{.78) 

398(.80) 

405(.81) 

375(.75) 

too  low 

344(.69) 

844(.69) 

54(.ll) 

76(.15) 

45(.09) 

104(.21) 

Table  25.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  UNMOD,  C=l/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.0872  .0864 

.0559  .0557 

.0701  .0683 

(.0049)  (.0048) 

(.0047)  (.0047) 

(.0238)  (.0213) 

3.0 

.0945  .0927 

.0754  .0745 

.0905  .0834 

(.0050)  (.0051) 

(.0033)  (.0035) 

(.0309)  (.0227) 

5.0 

.0855  .0819 

.0670  .0654 

.0851  .0736 

(.0106)  (.0112) 

(.0080)  (.0081) 

(.0369)  (.0230) 

7.0 

.0729  .0664 

.0524  .0503 

.0704  .0560 

(.0216)  (.0209) 

(.0105)  (.0106) 

(.0407)  (.0214) 

9.0 

.0578  .0502 

.0385  .0362 

.0551  .0397 

(.0317)  (.0285) 

(.0108)  (.0107) 

(.0430)  (.0184) 

11.0 

.0403  .0339 

.0273  .0251 

.0422  .0269 

(.0364)  (.0312) 

(.0098)  (.0095) 

(.0432)  (.0149) 

13.0 

.0290  .0240 

.0190  .0169 

.0322  .0178 

(.0355)  (.0298) 

(.0083)  (.0079) 

(.0430)  (.0115) 

15.0 

.0213  ,  .0175 

.0130  .0113 

.0247  .0116 

(.0330)  (.0275) 

(.0067)  (.0062) 

(.0427)  (.0085) 
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Table  26.  TWO-SIDED  90  %  COVERAGE  FRACTION(MOD  1,  C=  1/2) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

44(.09) 

30(.06) 

42(.08) 

26(.05) 

9(.02) 

8(.02) 

1.0 

cover 

444(.89) 

453(.91) 

445(.89) 

455(.91) 

451(.90) 

377(.75) 

too  low 

12(.02} 

17(.03) 

13(.03) 

19(.04) 

40(.08) 

115(.24) 

too  high 

55(.ll) 

39(.08) 

35(.07) 

29(.06) 

9{.02) 

7(.01) 

3.0 

cover 

439(.88) 

438(.88) 

449(.90) 

447(.89) 

457(.91) 

382(.76) 

too  low 

6(.01) 

23{.05) 

16(.03) 

24(.05) 

34(.07) 

111(.22) 

too  high 

100(.20) 

78(.16) 

39(.08) 

27(.05) 

9(.02) 

10(.02) 

5.0 

cover 

321(.64) 

340(.68) 

444(.89) 

448(.90) 

460(.92) 

354(.71) 

too  low 

79(.16) 

82(.16) 

17(.03) 

25(.05) 

31(.06) 

136{.27) 

too  high 

162(.32) 

134(.27) 

41(.08) 

26(.05) 

10(.02) 

10(.02) 

7.0 

cover 

194(.39) 

222(.44) 

443(.89) 

448(.90) 

456(.91) 

335(.67) 

too  low 

144(.29) 

144(.29) 

16(.03) 

26(.05) 

34(.07) 

155(.31) 

too  high 

224(,45) 

190(.38) 

43(.09) 

25(.05) 

10(.02) 

10(.02) 

9.0 

cover 

115(.23) 

149(.30) 

440(.88) 

440(.88) 

458(.92) 

324(.65) 

too  low 

161(.32) 

161(.32) 

17(.03) 

35(.07) 

32(.06) 

166(.33) 

too  high 

285(.57) 

217(.43) 

46(.09) 

24(.05) 

11(.02) 

10{.02) 

II.O 

cover 

51(.10) 

I19(.24) 

437(.87) 

435(.87) 

456(.91) 

307(.61) 

too  low 

164(.33) 

164(.33) 

17(.03) 

41(.08) 

33(.07) 

183(.37) 

too  high 

323(.65) 

250(.50) 

46(.09) 

19(.04) 

11(.02) 

11(.02) 

13.0 

cover 

13(.03) 

86(.17) 

437(.87) 

437(.87) 

454{.91) 

290(.58) 

too  low 

164(.33) 

164(.33) 

I7(.03) 

44(.09) 

35(.07) 

119(.40) 

too  high 

336(.67) 

290(.58) 

47(.09) 

16(.03) 

11(.02) 

12(.02) 

15.0 

cover 

O(.OO) 

46(.09) 

436(.87) 

434(.87) 

454(.91) 

274(.55) 

too  low 

164(.33) 

164(.33) 

17(.03) 

51(.10) 

35(.07) 

214(.43) 
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Table  27.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  MOD  1,  C=  1/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.1245  .1223 

.0959  .0949 

.2871  .2379 

(.0107)  (.0102) 

(.0080)  (.0077) 

(.1289)  (.1119) 

3.0 

.1783  .1647 

.1337  .1289 

.3048  .2088 

(.0341)  (.0249) 

(.0108)  (.0108) 

(.1359)  (.1049) 

5.0 

.1775  .1494 

.1251  .1164 

.2881  .1527 

(.0904)  (.0722) 

(.0223)  (.0223) 

(.1586)  (.1030) 

7.0 

.1535  .1267 

.1045  .0927 

.2657  .1040 

(.1083)  (.0864) 

(.0295)  (.0287) 

(.1775)  (.0902) 

9.0 

.1463  .1207 

.0830  .0695 

.2448  .0684 

(.1108)  (.0887) 

(.0320)  (.0299) 

(.1902)  (.0736) 

9.0 

.1452  .1197 

.0644  .0505 

.2276  .0442 

(.1113)  (.0891) 

(.0316)  (.0281) 

(.1991)  (.0578) 

9.0 

.1452  .1197 

.0495  .0361 

.2137  .0284 

(.1113)  (.0891) 

(.0295)  (.0249) 

(.2050)  (.0444) 

9.0 

.1452  .1197 

.0379  .0256 

.2023  .0182 

(.1113)  (.0891) 

(.0267)  (.0213) 

(.2089)  (.0336) 

57 


Table  28.  TW  O-SIDED  80  %  COVERAGE  FRACTION(MOD  1,  C=  1/2) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

68(.14) 

59(.l2) 

61(.12) 

54{.ll) 

28(.06) 

25(.05) 

1.0 

cover 

406(.81) 

396(.79) 

407(.81) 

405(.81) 

391(.78) 

322(.64) 

too  low 

26(.05) 

45(.09) 

32(.06) 

41(.08) 

81(.16) 

153(.31) 

too  high 

94(.19) 

75(.15) 

62(.12) 

47(.09) 

22(.04) 

20(.04) 

3.0 

cover 

381(.76) 

382(.76) 

404{.81) 

409(.82) 

406(.81) 

318(.64) 

too  low 

25(.05) 

43(.09) 

34(.07) 

44(.09) 

72(.14) 

162{.32) 

too  high 

129{.26) 

108{.22) 

64(.13) 

48(.10) 

20(.04) 

18(.04) 

5.0 

cover 

291(.58) 

302(.60) 

406(.81) 

405(.81) 

409(.82) 

305(.61) 

too  low 

80(.16) 

90(.I8) 

30(.06) 

47(.09) 

71(.14) 

177(.35) 

too  high 

19I(.38) 

164(.33) 

65(.13) 

48(.10) 

21(.04) 

17(.03) 

7.0 

cover 

165(.33) 

192(.38) 

406(.81) 

398(.80) 

401(.80) 

292(.58) 

too  low 

144(.29) 

144(.29) 

29(.06) 

54(.ll) 

78(.16) 

191(.38) 

too  high 

257(.51) 

217(.43) 

67{.13) 

47{.09) 

19(.04) 

18(.04) 

9.0 

cover 

82(.16) 

122(.24) 

405(.81) 

396(.79) 

403(.81) 

275(.55) 

too  low 

I61(.32) 

161(.32) 

28(.06) 

57(.ll) 

78(.16) 

207(.41) 

too  high 

307(.61) 

267(.53)  ; 

66(.13) 

45(.09) 

19(.04) 

17(.03) 

11.0 

cover 

29(.06) 

69(.14) 

405(.81) 

392(.78) 

403(.81) 

264{.53) 

too  low 

164(.33) 

164(.33) 

29(.06) 

63(.13) 

78(.16) 

219(.44) 

too  high 

334(.67) 

308(.62) 

69(.14) 

39(.08) 

19{.04) 

13(.03) 

13.0 

cover 

2(.00) 

28(.06) 

402(.80) 

396(.79) 

401(.80) 

253(.51) 

too  low 

164(.33) 

164(.33)  ! 

29(.06) 

65(.I3) 

80(.I6) 

234(.47) 

too  high 

336(.67) 

333(.67)  1 

70(.14) 

34{.07) 

19(.04) 

12(.02) 

15.0 

cover  j 

O(.OO) 

3  (.01) 

400(.80) 

400(.80) 

403{.81) 

245(.49) 

too  low 

1 

164(.33) 

164(.33)  i 

30(.06) 

66(.13) 

78(.16) 

243{.49) 

Table  29.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  MOD  1,  C=l/2) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

1.0 

.0968 

.0957 

.0746 

.0742 

.2321 

.1925 

(.0082)  (.0080) 

(.0062)  (.0061) 

(.1207)  (.0972) 

3.0 

.1359 

.1296 

.1032 

.1010 

.2440 

.1648 

(.0240) 

(.0200) 

(.0085)  (.0085) 

(.1349)  (.0879) 

5.0 

.1311 

.1181 

.0953 

.0912 

.2284 

.1193 

(.0655) 

(.0574) 

(.0175)  (.0175) 
_ 1 

(.1592)  (.0843) 

7.0 

.1126 

.1002 

.0781 

.0726 

.2069 

.0807 

(.0784) 

(.0686) 

(.0229)  (.0226) 

(.1748)  (.07.26) 

9.0 

.1073 

.0955 

.0607 

.0545 

.1865 

.0528 

(.0804)  (.0704) 

(.0245)  (.0236) 

(.1829)  (.0586) 

11.0 

.1065 

.0947 

.0459 

.0396 

.1697 

.0340 

(.0808)  (.0707) 

(.0237)  (.0221) 

(.1876)  (.0456) 

13.0 

.1065 

.0947 

.0343 

.0283 

.1566 

.0218 

(.0808) 

(.0707) 

(.0217) 

(.0196) 

(.1906) 

(.0347) 

15.0 

.1065 

.0947 

.0256 

.0200 

.1463 

.0139 

(.0808)  (.0707) 

(.0192)  (.0167) 

(.1924)  (.0261) 

59 


Table  30.  TWO-SIDED  90  %  COVERAGE  FRACTION(MOD  1,  C=  1/10) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

34{.07) 

26(.05) 

31(.06) 

24(.05) 

I4(.03) 

10(.02) 

1.0 

cover 

453(.91) 

454(.91) 

448(.90) 

443(.89) 

454(.91) 

440(.88) 

too  low 

13(.03) 

20(.04) 

21(.04) 

33(.07) 

32(.06) 

50(.10) 

too  high 

24(.05) 

16(.03) 

27(.05) 

21(.04) 

20(.04) 

13{.03) 

3.0 

cover 

460(.92) 

456(.91) 

450(.90) 

1  444(.89) 

463(.93) 

448{.90) 

too  low 

16{.03) 

28{.06) 

23(.05) 

35(.67) 

I7(.03) 

39(.08) 

too  high 

37(.07) 

20(.04) 

27(.05) 

19(.04) 

19{.04) 

11(.02) 

5.0 

cover 

453{.91) 

452(.90) 

452(.90) 

445(.89) 

458{.92) 

446(.89) 

too  low 

10(.02) 

28(.06) 

21(.04) 

36(.07) 

23(.05) 

43(.09) 

too  high 

37(.07) 

25(.05) 

27(.05) 

18(.04) 

2I(.04) 

10{.02) 

7.0 

cover 

438(.88) 

434(.87) 

454(.91) 

442{.88) 

455(.91) 

439(.88) 

too  low 

25(.05) 

41{.08) 

19(.04) 

40{.08) 

24(.05) 

51(.10) 

too  high 

47(.09) 

27(.05) 

26(.05) 

16(.03) 

23{.05) 

8(.02) 

9.0 

cover 

357(.71) 

377(.75) 

455(.91) 

438(.88) 

454(.91) 

435(.87) 

too  low 

96(.19) 

96(.19) 

19(.04) 

46(.09) 

23(.05) 

57(.ll) 

too  high 

59(.12) 

29(.06) 

26(.05) 

12(.02) 

23(.05) 

6(.01) 

11.0 

cover 

233(.47) 

263(.53) 

455(.91) 

434(.87) 

453(.91) 

424(.85) 

too  low 

208(.42) 

208{.42) 

19(.04) 

54(.ll) 

24(.05) 

70(.14) 

too  high 

59(.12) 

38(.08)  1 

26<.05) 

n{.02) 

23{.05) 

5(.01) 

13.0 

cover 

154(.31) 

175(.35) 

454{.91) 

433(.87) 

453(.91) 

419(.84) 

! 

too  low 

287(.57) 

287(.57) 

20(.04) 

56(.ll) 

24(.05) 

76(.15) 

too  high 

83(.17) 

35(.07) 

26(.05) 

7(.01) 

23(.05) 

3(.01) 

15.0 

cover 

73(.15) 

121(.24) 

456{.91) 

436{.87) 

453(.91) 

418(.84) 

too  low 

344(.89) 

344(.89) 

18(.04) 

57(.ll) 

24{.05) 

79(.16) 
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Table  31.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(90% 
C.I,  MOD  1,  C=l/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.1121  .1105 

.0718  .0714 

.1276  .1187 

(.0064)  (.0062) 

(.0061)  (.0060) 

(.0812)  (.0704) 

3.0 

.1223  .1185 

.0972  .0953 

.1431  .1222 

(.0062)  (.0066) 

(.0042)  (.0044) 

(.0762)  (.0457) 

5.0 

.1123  .1046 

.0871  .0837 

.1300  .1021 

(.0132)  (.0142) 

(.0101)  (.0103) 

(.0765)  (.0379) 

7.0 

.0988  .0847 

.0688  .0644 

.1077  .0761 

(.0285)  (.0267) 

(.0134)  (.0135) 

(.0760)  (.0325) 

9.0 

.0806  .0640 

.0512  .0464 

.0861  .0533 

(.0437)  (.0363) 

(.0139)  (.0137) 

(.0749)  (.0268) 

11.0 

.0573  .0432 

.0369  .0321 

.0681  .0360 

(.0513)  (.0398) 

(.0128)  (.0122) 

(.0733)  (.0211) 

13.0 

.0416  .0305 

.0261  .0217 

.0540  .0237 

(.0506)  (.0380) 

(.0109)  (.0101) 

(.0716)  (.0160) 

15.0 

.0306  .0223 

.0182  .0145 

.0433  .0154 

(.0472)  (.0350) 

(.0089)  (.0079) 

(.0701)  (.0118) 

) 
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Table  32.  TWO-SIDED  80  %  COVERAGE  FRACT10N(M0D  1,  C=  1/10) 


TIME 

cov- 

K.M.E 

M.L.E 

A.R.E 

ERAGE 

LOG 

ASIN 

LOG 

ASIN 

LOG 

ASIN 

too  high 

61{.12) 

47(.09) 

56(.ll) 

48(.10) 

35(.07) 

28{.06) 

1.0 

cover 

404(.81) 

402(.80) 

400(.80) 

403(.81) 

396(.79) 

392(.78) 

too  low 

35(.07) 

51(.10) 

44(.09) 

49(.10) 

69{.14) 

80(.16) 

too  high 

49(.10) 

38(.08) 

52(.10) 

42(.08) 

49(.10) 

37(.07) 

3.0 

cover 

415(.83) 

417(.83) 

401(.80) 

403(.81) 

408(.82) 

399(.80) 

too  low 

36(.07) 

45(.09) 

47(.09) 

55(.ll) 

43{.09) 

64{.13) 

too  high 

56(.ll) 

47(.09) 

52(.10) 

40(.08) 

49(.10) 

39(.08) 

5.0 

cover 

411(.82) 

392(.78) 

397(.79) 

395(.79) 

409(.82) 

390(.78) 

too  low 

33(.07) 

61(.12) 

51(.10) 

65(.13) 

42(.08) 

71(.14) 

too  high 

55(.ll) 

43(.09) 

51(.10) 

39(.08) 

47{.09) 

34{.07) 

7.0 

cover 

410(.82) 

390(.78) 

395(.79) 

396(.79) 

407{.81) 

392(.78) 

too  low 

35(.07) 

67(.13) 

54(.ll) 

65(.13) 

46(.09) 

74(.15) 

too  high 

71(.14) 

53(.ll) 

51(.10) 

34(.07) 

48(.10) 

31(.06) 

9.0 

cover 

333(.67) 

351(.70) 

394(.79) 

399(.80) 

407(.8I) 

384(.77) 

too  low 

96(.19) 

96(.19) 

55(.ll) 

67(.13) 

45(.09) 

85(.17) 

loo  high 

79(.16) 

59(.12) 

53(.ll) 

31(.06) 

47(.09) 

25(.05) 

11.0 

cover 

213(.43) 

233(.47) 

392{.78) 

400(.80) 

407(.81) 

381(.76) 

too  low 

208(.42) 

208(.42)  1 

55(.ll) 

69(.14) 

46(.09) 

94(.19) 

too  high 

84(.17) 

56(.ll) 

55(.ll) 

29(.06) 

47(.09) 

24(.05) 

13.0 

cover 

129(.26) 

157(.31)  i 

391(.78) 

398(.80) 

406(.81) 

379(.76) 

too  low 

287(.57) 

287(.57) 

54(.ll) 

73(.15) 

47(.09) 

97C.I9) 

too  high 

110(.22) 

58(.12) 

54(.ll) 

26{.05) 

48(.10) 

21(.04) 

15.0 

cover 

46(.09) 

98(.20) 

392(.78) 

398(.80) 

405(.81) 

373(.75) 

too  low 

344(.69) 

344(.69) 

54(.ll) 

76(.15) 

47(.09) 

106(.21) 
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Table  33.  AVERAGE  AND  STANDARD  DEVIATION  OF  HALF  LENGTH(80% 
C.I,  MOD  1,  C=  1/10) 


TIME 

K.M.E 

M.L.E 

A.R.E 

LOG  ASIN 

LOG  ASIN 

LOG  ASIN 

1.0 

.0872  .0864 

.0559  .0557 

.0981  .0908 

(.0049)  (.0048) 

(.0047)  (.0047) 

(.0665)  (.0563) 

3.0 

.0945  .0927 

.0754  .0745 

.1070  .0929 

(.0050)  (.0051) 

(.0033)  (.0035) 

(.0593)  (.0357) 

5.0 

.0855  .0819 

.0670  .0654 

.0952  .0776 

(.0106)  (.0112) 

(.0080)  (.0081) 

(.0574)  (.0293) 

7.0 

.0729  .0664 

.0524  .0503 

.0771  .0578 

(.0216)  (.0209) 

(.0105)  (.0106) 

(.0571)  (.0250) 

9.0 

.0578  .0502 

.0385  .0362 

.0600  .0404 

(.0317)  (.0285) 

(.0108)  (.0107) 

(.0561)  (.0205) 

11.0 

.0403  .0339 

.0273  .0251 

.0461  .0273 

(.0364)  (.0312) 

(.0098)  (.0095) 

(.0548)  (.0161) 

13.0 

.0290  .0240 

.0190  .0169 

.0353  .0180 

(.0355)  (.0298) 

(.0083)  (.0079) 

(.0535)  (.0122) 

15.0 

.0213  .0175 

.0130  .0113 

.0274  .0116 

(.0330).  (.0275) 

(.0067)  (.0062) 

(.0524)  (.0089) 
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